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1 Introduction 



Invariant group integrals are encountered in many problems in physics. Examples are the 
one-link integrals in mean field or strong coupling expansions in Euclidian lattice gauge the- 
ories [T] or mass-gap calculations in the Hamiltonian formulation of these theories [2j. They 
are used in the exact solution of two-dimensional lattice gauge theories [3], the matching 
of gauge theories to chiral models [4] and the loop formulation of quantum gravity on spin 
network states [5]; they appear in random matrix theory [6j and its widespread applications 
in nuclear physics [7], quantum chaos and transport in mesoscopic devices |8j as well as 
in quantum information theory [9]. Apart from that, they are tightly connected to vari- 
ous enumerative problems in mathematics such as counting the number of invariants in a 
given tensor product of group representations or the number of Young tableaux of bounded 
height [lO]. 

We consider invariant integrals over compact Lie groups with left-right invariant Haar mea- 
sure (i/iHaar- For example, integrals of the form 

ZgUJ^) = J df^Us^Ua) exp (tr (j^g) + tr (jg^)) (1) 

with a matrix valued source j can be calculated for the U{N) groups [H]. The solution 
of the Schwinger-Dyson equations corresponding to the left and right group actions on the 
source yields the closed form expression in terms of Bessel functions, 

where {zb/2)'^ are the left- and right-invariant eigenvalues of j ■ j^. For the SU{N) groups 
there are further invariants independent of j\ for example det j, and no comparable simple 
solution to the Schwinger-Dyson equations are known. 

An alternative approach based on an explicit parametrization of the group elements has 
been proposed in [12]. Here one ends up with series representations of the form 

Zg{J,J^)= Y1 (3) 

0<ni,...,ns<oo 

where the algebraically independent invariants, Xk{gjg') = Xk{j). This 

method is only applicable to small groups with sufficiently simple parametriziations. For 
example, for SU (2) one obtains the series 

oo _ 
n=0 ^' 
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Already for SO (3), the result is a triple series expansion in powers of the three independent 
invariants. For SU{3), a convenient parametrization of the group elements leads to the 
series [T2l 



4 

Zsui3)ij,J^)= (ni + 2n2 + 3n3 + n4 + 2) ! (722 + 2^3 + ^4 + 1) ! n 

ni,...,n4 ^ ' ^ ' p=l 

with four left-right invariants 

xi = tr (jj^), a;2 = I (a;i - tr (jj^y) , = det{jj^), = detj + detj^. (6) 

Since X4 is not invariant under U{3) transformations one finds a triple series for this group. 
The corresponding results for U{2) and U{3) actually coincide with the expression in terms 
of Bessel functions, eq. ([2|) [12]. 

Numerous results for so-called n-vector integrals containing only n columns of unitary ma- 
trices were derived in [13] and more recently in [H] using an elegant method based solely on 
the unitary constraint and left-right invariance of the Haar measure. These methods were 
applied to orthogonal groups in |15J. 

In ([H), is a generating function for invariant integrals of arbitrary functions on the group. 
In many applications one only needs integrals of class functions. Such functions are constant 
on conjugacy classes, 

F{ghg-') = F{h), (7) 

such that we can consider them as functions of the maximal Abelian torus in G. By the 
Peter- Weyl theorem, the group characters form an orthonormal basis on the Hilbert space 
of square integrable class functions. Characters can be computed with the help of Weyl's 
character formula; they are polynomials of the 'fundamental characters' Xp belonging to the 
fundamental representations with highest weights p = 1, . . . , r. 

One of the authors was involved in the mean field analysis of effective models for pure gauge 
theories near their phase transition point [16]. There one is confronted with calculating 
invariant integrals of the type 

Zg{u) = / dfi,cd{g) exp {^UpXpig)), (8) 

p=i 

where dfi^ed is the reduced Haar measure on G. Such invariant integrals together with their 
Ward identities also show up in inverse Monte-Carlo simulations where one calculates the 
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couplings in the Polyakov loop dynamics of gauge theories [T7]. Other applications are, e.g., 
the computation of glueball masses in Hamiltonian SU (N) lattice gauge theories [21 [18] 
or studies of the strong coupling limit of these theories |T9]. The function Zg{u) is the 
generating function for the moments of all fundamental characters, 

ami+...+mr r 

for r-tuples m = (ui, . . . , 

In the literature there seem to be no suitable Ward identities for Zg{u) for arbitrary 
groups G. Thus we decided to publish our findings since they could be useful for colleagues 
confronted with similar invariant integrals. 

The paper is organized as follows: The next section gives an overview of the known results 
for the generating functions Zg{u) for G = U{N) and G = SU{N) and extend those 
for SU{N) in the special case of a generating function for the defining representation. In 
section [3] we derive what we call geometric Ward identities, since they are based on the 
invariance of the Haar measure. These results hold for all compact and simple Lie groups. 
In section [4] our method is applied to the simple rank 2 groups SU{3), Spin(5) and G2 and 
in section [5] to 5't/(4). In section [6] an alternative and more analytic method is presented 
which sheds further light on the properties of the reduced Haar measures. It turns out that 
the square of the Jacobian of the transformation ^ ^ X from the angles parametrizing the 
maximal Abelian torus to the fundamental characters is proportional to the density of the 
reduced Haar measure. We have no general proof of this conjecture but have checked it for 
groups with rank 2 and 3. Based on the conjecture we find alternative Ward identities which 
are applied to the groups with rank 2. In section [7] we use both types of Ward identities 
to derive recursion relations for the moments tmi,...,mr in ®- The appendices contain a 
detailed description of the solution of the Ward identities derived in section 14.11 as well as 
tables of the lowest moments for the above rank 2 groups. In the conclusions we comment 
on possible applications of our Ward identities. 



^Note that for j — uil, eq. (HJ yields the generating function for the defining representation of U{N) 
and its complex conjugate; however, for iV > 3 there is no direct way to derive the generating function 
^u{N){u) for all fundamental characters from the generating function Zij(^j^-^{j,j''). 
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2 Results for U{N) and SU{N) 

Let z be an arbitrary element of the group center. By Schur's lemma it acts on Xp(fl') by 
multiplication with a factor Zp such that for the r-vector x = (Xi? • • • ? Xr)*, 



= with D{z) = diag{z,,...,Zr). (10) 

By invariance of the Haar measure under g i-^ zg, this implies the symmetry 

Zg{u) = Zg {D-\z)u) G center. (11) 

This observation proves to be crucial for explicit computations of Zq in the following sec- 
tions. We now briefly summarize the known results for G = U{N) and extend those for 
SU (N) in the special case of a generating function for the defining representation. 
For U{N) the fact that the reduced Haar measure factorizes into a flat measure times the 
absolute square of a Vandermonde determinant detA{{p), 

dfiredig) = det A{ip) det A*{if) d'^p, (12) 

facilitates the integration such that the generating function 

jd,,^eMnXrg + virg') (13) 



'U{N) 



for integrals over the characters of the deflning representation and its complex conjugate 
can be computed explicitly. An integration over the maximal torus parameterized by the 
angular variables ipi, . . . ,{pj^ yields the closed expression 



^u{N)iu,v) = det 



/ /o h ■■■ In-i\ 
h Iq ■ ■ ■ 

\ In-1 In-2 ■ ■ ■ Iq / 



{2^/mJ) . (14) 



In accordance with the general result (flTl) . this function is invariant under center transfor- 
mations, 

ZuiN) (e'"M,e-*"^;) = ZuiN){u,v). (15) 

It can easily be checked that for v = u* this function is the limit of ([2]) for Zf, 4m*m; this 
limit has to be taken with caution since both the numerator and denominator in (12]) vanish 



4 



for two or more coinciding eigenvalues such that I'Hopital's rule and certain Bessel function 
identities are needed. 

In order to calculate the analogous generating function for SU{N) we follow [20] and insert 
the constraint det g = 1 on the maximal Abelian torus in the form 



•+Viv) 



into the invariant integral (|T3|) and find 

/ In 



Zsu{N) (m, f ) = 2^ I - I det 

ngZ 



In+1 

L 



In+N-l \ 
In+N-2 

In J 



(16) 



\ In-N+1 In~N+2 ' ' ' 1 

see [IH]. The generating function for SU{N) is left invariant by center transformations 
for which 



UV — i> uv, 



u 

V 



V 



In the case of S'f/(2), the sum can be worked out and leads to 

Zsu(2){u,v) = qFi[2 \w'^] = /o(4w) - hi^w), w 



U + V 



(17) 



(18) 



The case of SU{3) is the subject of section [4Tl We are not aware of similar explicit results 
for SU{N) with > 4. In these cases, the complexity increases since the complex conjugate 
fundamental representations are inequivalent. 

In (fTSl) . the hypergeometric function o-^i does not appear accidentally; for all SU{N) the 
generating function for the defining representation, 

(19) 



Zsu(N)iu) = J dfij-edig) exp{utvg) 
is one of the generalized hypergeometric functions 



p-Tq 



tlx , . . . , dp 



X 



Ex 



n=0 



an 



in + ai 



in + 



in + ha 



(20) 



with ao = 1. In the course of the paper, we will mostly identify them as solutions to the 
generalized hypergeometric differential equation, 



q p 

ell{9 + k-i)~xl[{ 

i=l 1=1 



p-Tq 



ai. 



, Clp 
,bq 



X 



with 6 = x^. (21) 

ax 



For g = (or p = 0), the first (or second) products of differential operators are to be 
replaced by the identity operator. 

Center symmetry entails that the function Zsu{n){u) in (fTOll is in fact only a function of 
X = det(Ml) = , and Zsu{n){u) = Z{x) satisfies the differential equation 

^ (x^-'Zix)) = Z{x). (22) 

The solution is the hypergeometric function qFn-i [2, 3, ... , A^|x] such that 

ZsuiN)iu) = oFN^i[2,3,...,N\u^]. (23) 

Since for SU{2)^ tig = ti g\ this generalizes the standard result (fTHl) (with u = v) hy 
Arisue [2T]. Eq. (l23l) follows from (flGl) when v tends to zero so that 

ZsuiN){u) = J2^''-detA^-\ (A("))^=|^^ forn + g-p>0, _ ^^^^ 
n>o I else 

If we multiply the p'th row of A*^") with (n + A^ — p)!, it is easy to calculate the determinants 
of these Toeplitz matrices, 

dot A^"^ - TT ^ =^ = (2^) 

A^-L(n + p)! detAW {72 + 1) ■ ■ ■ {n + N)' ^ ' 

This proves eq. (l23l) . 



3 Ward identities for generating functions 

We denote the left derivative in the direction of the Lie algebra element by La, i.e. 
Laf{g) = ^|(=o/( exp(itT(j)5f) for some function / on G. The Haar-measure is left (and 
right) invariant, thus 

j dfinUg) {LJ){g) = 0, / G L^{G). (26) 
For class functions F and F the function 

^La[F- LaF^ = L [fLF^ = FL^F + LF ■ LF (27) 

a 

is a class function as well. In order to see this, we may assume that F and F are basis 
elements (i.e., characters x^i and of some representations with highest weights /i, u). In 
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J2a(Xf^ ' ^aXu + LaX^i ' LaXu), the first part of the sum is a class function since the quadratic 
Casimir operator Ylia'^a commutes with all group elements, and the second part is a class 
function since invariance of the Killing metric tr TaT^, under adjoint action by some group 
element h implies that hTah~^ can be expanded as Ra^T^ with an orthogonal matrix R. 
Thus, (l26ll with f = F ■ LaF reduces to an integral over the maximal Abelian torus, 







J dfi,^d L (^FLF^ = j dfired [fL^F + LF ■ LF^ . (28) 



We take F to be a fundamental character Xp with p G {1, . . . , r}. The xi? ■ ■ ■ ^Xr are good 
coordinates for the maximal torus such that any class function can be thought of as a 
function of these characters, F = F{xiy • • • ? Xr)- Then the identity (l28ll reads 

= I t//i.ed (^Fix)L\, + J2(^Lxp)^^iLx,?j . (29) 

Every character x^l of a representation with highest weight /i is an eigenfunction of the 
Laplacian with eigenvalue — c^, where is the value of the quadratic Casimir in V^, 

L^Xf^ = -C/.Xm- (30) 

To calculate the last term in (l29ll we decompose the tensor product of Vu®V,, into irreducible 
piecesp 

^M®K = 0C^.^A, such that x,^Xu = J2C^,Xx- (31) 

A 

Acting with on this relation and using (l30l) we find the useful relation 

{LXf,) ■ {Lxu) = ^(c^ + c^)Xf.Xu - ^ 5Z C-^i. cxxx- (32) 

A 

with Clebsch-Gordan coefficients C^^ and second order Casimirs Cn- Now we may rewrite 



the Ward identity ( l29ll as follows, 



Cp J dflredXpF{x) = ^5^(Cp + Cg) j df^redXpXg 



dF{x) 

dXg 



\ Yl ^ii / ^-"-d XA ^1^, 1 < p < r. (33) 



^The characters x\ are polynomials x\ — Xa(xi, • • ■ , Xr) of the fundamental characters. For groups of 
higher rank it can be cumbersome to calculate these polynomials. 
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We choose the class function F = exp{u ■ x) such that 

f dF dF 

J dfiredF = Zg{u), — = UqF and — = XgF. (34) 
Then, the identity fl33ll translates into the following master equation for Zq, 



-^^CpqCxUqX\{d)ZG{u), p=l,...,r, (35) 

q,X 

where Xa(^) is the differential operator obtained by formally evaluating the polynomial 
X\{x) cit {d/dui, . . . , d/dur). The properties of the group G enter these Ward identities at 
three places: via the polynomials Xx{x)i the values of the quadratic Casimir operators 
and the Clebsch-Gordan coefficients Cp^. The coefficient functions of these linear partial 
differential equations are constant or linear functions of the variables ui,...,Ur. Their 
complexity depends crucially on the polynomials xa- For SU{3) and Spin(5) all xx 
quadratic, and one ends up with second order differential equations. 

For explicit calculations, we have to fix our Lie algebra conventions which were chosen 
to allow for an easy comparison with the computer algebra program LiE [22]. In these 
conventions, the Cartan matrix is given by 

2(a(p),a(g)) 

^pq = 7 T (36) 

("(<?), "(g)) 

in terms of the simple roots . . . , a(r). Thus, simple roots and fundamental weights are 
connected by the relation 

q 

The shortest simple root has squared length 2. In particular, for simply laced groups, 
the Cartan matrix reduces to Kpg = (a(p),a(q)). Arbitrary roots and weights are linear 
combinations of the simple roots and fundamental weights, respectively, 

a = mpa(p) = [mi, . . . , irir] and /i = ^pAt(p) = [ni, . . . , n^]. (38) 

The Weyl vector 

1 

c(>0 p=l 
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plays an important role in the theory of representations; we will mostly need it for its 
appearance in the formula for the value of the second order Casimir operator in a represen- 
tation V^, 

= (/i,/x + 2p). (40) 

With the help of ( i37l ) and = 6pqVp the Casimir of the representation with highest 

weight /i = [rii, . . . ,nr] can be written as 

= (n, K~^n') with n'^ = Vp{np + 2). (41) 

For simply laced groups all Vp are equal to 1 so that = + 2. 

4 Ward identities for groups of rank 2 

In this section, we are going to investigate and exploit the geometric Ward identities (i35l l 
for the groups SU{3), Spin(5) and G2- For these rank 2 groups, the generating function ([8]) 
depends on two variables ui = u and U2 = v. 

4.1 The group SU{3) 

As a simply-laced group SU (3) has a symmetric Cartan matrix, 

(42) 

and the quadratic Casimir of the representation with highest weight /i = [^1,^2] is 

2 

C/, = - {nj + nl + ^1^2 + 3ni + 3n2) . (43) 

The fundamental 3-dimensional representation 3 = [1, 0] and its complex conjugate 3 = [0, 1] 
both have Casimir 8/3. Since 

3®3 = 6©3, 3(g)3 = 6©3, 303 = 1© 8, (44) 

the A-sum in (l35l) contains both fundamental, two sextet and the octet representations. 
The singlet representation has vanishing Casimir invariant and does not contribute. The 
Casimir operator on the sextets 6 = [2, 0] and 6 = [0, 2] takes the value 20/3, and the octet 
8 = [1,1] has Casimir 6. To derive explicit Ward identities we must express the characters 
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of the representations 6, 6 and 8 in terms of the fundamental characters. By (ISTll . eq. (I44l) 
yields 



X6 = xl- Xs, X6 = xl- Xs' X8 = X3X3 - 1- (45) 
Thus the differential operators X\{x) iii the Ward identity ( l35l) read 

X3 = du, X3 = dv: Xe = dl- d^, XG = dl- <9„, xs = dud^ - 1, (46) 
and the two Ward identities in (l35ll for the generating function 

Zsu{3){u,v) = dfired{g)ew{uX3{g) + VX3{9)) with x3{g)=tT^g (47) 



have the simple for: 




{2udl + vdud, + 8du-6ud,-9v)Zsui3) = 0, (48) 
{2vd^^+ud,du + 8d,-6vd^-9u)Zsui3) = 0. (49) 

Note that by the property dfi{g) = dfi{g^^) of the Haar measure, Zsu{3){u,v) = Zsu{3)iv,u) 
is a symmetric function. Thus it is sufficient to study one of the two identities. In the 
defining representation of S'f/(3), the center Z3 acts by multiplication with exp(27r2/3)l, 
and the symmetry (fTT]l reads 

Zsu(3) (e'-'/^M, e-^^^'S) = ZsuiS){u, v). (50) 

Together with the symmetry in its two arguments this suggests that Zsu(3) is just a function 
of the combinations v? + f ^ and uv. In fact, we prove in the appendix that the solution to 
(Il8l) is given by 

We have found only one comparably simple series representation for the generating function 
in the literature [18]. To arrive at their result the authors took the quartic series ([5]) and 
calculated 2 of the 4 infinite sums. It seems evident that our method based on geometric 
Ward identities is more efficient to find simple series representation for generating functions. 

^The action of left-derivatives on characters can also be worked out concretely in a matrix representation 
although this requires more computational effort, cf. }23j . 
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As two special cases, we will restrict Zsu{3){u,v) to the w-axis (i.e., v = 0) and the diagonal 
(i.e., u = v). Instead of performing a resummation of the general result (ISTl) . we derive 



differential equations for Zsu{3) in these cases. First we solve the Ward identities (148|49p on 
the M-axis, where 

{2udl + 8duZ -Qud^)Zsui3)\v=o = 0, {ud^du + 8d^ - 9u)Zsu{3)\v=o = 0- (52) 



To get rid of the f-derivatives at f = we act with udu on the first equation and use (l52l) 
to eliminate the terms d^Zsu{3) and dudyZsu{3)- We find the ordinary differential equation 

{u^dl + Uudl + 28du - 27u^) Zsu{3) («, 0) = 0. (53) 

Since Zsu(3)iu,Q) only depends on v? we may equally well use x = as a new variable, 
Zsu{3){u, 0) = Ysu{3){x). Then the differential equation takes the simpler form 

{x^dl + &xdl + 69, - 1) Ysui,3) = {x^Ysui:3))"' - Ysu(3) = 0, (54) 

which has a solution in terms of a generalized hypergeometric function (cf. (l22l) ). Ysu{3){,x) = 
0-F2 [2, 3 |x] . This is just the result (l23l) for the group SU{3). 

In order to find a differential equation on the diagonal u = v, we act by the operator 
{3udu + 5ud^ + 5 + 4u) on (Il8l) and evaluate the result at u = v =: t. In terms of 

d (f 

-^Zsu{3)it,t) = {du + dy)Zsu{3), -^Zsu{3)it,t) = C^duu + 2d^y)Zsu{3), 



'j;^Zsu{3)it,t) = {2duuu + Qduuv)Zsu{3), (55) 



the resulting ordinary differential equation reads 

d^ , d"^ . c? 



{f— + t(10 - t)— - 2{l2e + t - 10)- - m(3t + h))Zsui3) = 0. (56) 
This is solved by a function 

00 

Zsui3){t.t) = l + e + Y,'^t'' (57) 

n=3 

where the coefficients satisfy the recursion relation 

On+i = 7 — — nr? — TTT (^(^ + l)"^" ~ I2n{2n + 3)a„_i - 36n(n - l)a„_2) = 0. (58) 
[n + A)[n + b) 

Together with oq = 1, ai = 0, and 02 = 2, this determines all coefficients in the expan- 
sion (1581). 
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4.2 The group Spin(5) 

Since Spin(5) is not simply-laced, the Cartan matrix 

^spi„(5) = ( ) (59) 

is not symmetric. With the convention in (l36l) . we take ai to be the longer root, = 4 and 
= 2. The eigenvalue of the Casimir operator of the representation with highest weight 

fj, = [rii, n2\ reads 

= 2n\ + n\ + 2nin2 + 6ni + An2. (60) 

The fundamental representations are the 5*0(5) vector representation 5 = [1,0] and the spin 
representation 4 = [0, 1]. The center Z2 is generated by —1 in the spin representation and 
acts trivially in the vector representation; the center symmetry (fTTl) implies that 



^Spin(5)(M, ^) = y dUredig) GXp {uX^ig) + VX^ig)) = ^Spm(5) (^i, -v) (61) 

is an even function in v. For the geometric Ward identities we need the tensor products 

5® 5 = 10 100 14, 4®4=1©5©10, 5 ® 4 = 4 © 16. (62) 
With 10 = [0, 2], 14 = [2, 0] and 16 = [1, 1], we obtain for the Casimir operators (liTD : 

C5 = 8, C4 = 5, cio = 12, Ci4 = 20, C16 = 15. (63) 
Together with (l62ll . (l3Tl) implies that 

Xio = X4 - X5 - 1, Xi4 = xl-xl + X5, X16 = X5X4 - X4, (64) 
which leads to the Ward identities 

{uiAdl-2dl-Ad^ + Q)+vi5d,-dM-8d^)Zspin{5)iu,v) = 0, (65) 
{u{5dy - dudv) + v{2du - + Q) - 59„) ^spm(5) (u, v) = 0. (66) 

Since the center of Spin(5) is smaller than the center of SU (3) these differential equations are 
more complicated than the corresponding SU{3) equations in (Il8l) and fllQll . The character- 
istics of the second (generally hyperbolic) equation are given by m = const, and ^ = const., 
respectively. These families of lines coincide for m = (where ( l66l) is parabolic), and we 
may solve the characteristic problem given a solution of both equations for m = 0. 
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The restriction of the Ward identities to the f -axis 

(89„ - v{hd, - 9„5,))Zspin(5)L=o = 0' (5^- - - dl + 6))Zspi„(5)L^o = 0, (67) 

can be solved if we differentiate the second equation with respect to v and use the two 
relations (l67D to eliminate the w-derivatives at w = 0. We find the following ordinary 
differential equation for Zspin(5)(0, f ), 

{v'^dl + l?>vdl - IQv^d, + 359, - 48i;) Zspin(5)(0, v) = 0. (68) 

Since 2'spin(5) is an even function in v we set -Z'spm(5)(0, f ) = Zq{x = f^) and obtain the 
simple equation 

[x^dl + Sxdl - Axd^ + 129^ - 6) Zq{x) = 0. (69) 
The solution is a hypergeometric function 1-^2 [3^4 |4a;] so that 

^spin(5)(0, v) = j dfx,,d{g) exp (vxiig)) = 1^2 [f} \^v^] . (70) 
Plugging this into the second Ward identity ( l66l ) we obtain a recursive solution 

Zs,H^)iu,v) = f^^^^^^^u-, (71) 

n=0 

where 

Zi(x) = {6d^ + 2xdl-3)Zo{x), 

Z^{x) = [(5 + n)9, + 2a;a2-3]Z,_i(x)-5(n-l)9,Z„_2(x) for n>2. (72) 
As a special case, let us consider the Spin(5)-Ward identities on the w-axis, 

{u{2dl -dl- 2du + 3) - 49„)Zspin(5)L=o = 0' ("(59, - dM - 59,)Zspin(5)L=o = 0' C^S) 
We differentiate the first equation with respect to u and obtain 

{2dl - hdl - 29„ + 3 + u{2dudl -dl- 2dl + 39„)) Zspi„(5) = 0. (74) 

Eq. (1731) together with the v-derivative of (l66l) at w = can be used to eliminate the 
w-derivatives, 

{u^dl + (10m - ?>u^)dl + (20 - 12n - l?>u^)du - 30u + iW) Zspin(5) (m, 0) = 0. (75) 
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Very probably this cannot be converted into a differential equation for a generalized hyper- 
geometric series. The differential equation is solved by the power series 



oo , 

h 



^Spin(5) (m, 0) = 1 + ^2 + ^ ^ (76) 



n=3 



provided that the coefficients satisfy the recursion relation 

bn+i= , ,Z ^ ^J 3n{n + 3)bn + n{13n + 17)bn-i-15n{n-l)bn-2), n>2. (77) 
Together with bo = 1, bi = and 62 = 2, this determines all coefficients in the expansion 

m. 



4.3 The group G2 

For this group, the Cartan matrix is 

^°'K-3 "2')- 

Hence, ai is a short root and 02 a long root, af = 2 and a"^ = 6. The quadratic Casimir of 
the representation with highest weight fi = [^21,^2] reads 

Cf, = 2nl + 6nl + GniKa + lOni + 18n2. (79) 

The first fundamental representation 7 = [1,0] coincides with the subrepresentation of the 
8-dimensional spinor representation of Spin(7) leaving an arbitrary spinor fixed, and the 
second is just the adjoint representation 14 = [0,1]. The center of G2 is trivial and we 
expect no symmetries of the generating function. 
For the Ward identity we need the tensor products 

7® 7 = 1© 7© 14027, 

7 ©14 = 7 ©27 ©64, 

14 ©14 = 1© 14 ©27 ©77 ©77', 

7 © 27 = 7 © 14 © 27 © 64 © 77'. (80) 

We need the last product in order to express the characters xx as functions of the funda- 
mental characters. Note that two irreducible representations of dimension 77 appear in the 
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decompositions. We identify 27 = [2,0], 64 = [1, 1], 77 = [0,2] and 77' = [3,0] so that the 
quadratic Casimir operators (|4T1) take the following values: 

C7 = 12, ci4 = 24, C27 = 28, C64 = 42, C77 = 60, C77/ = 48. (81) 

Furthermore, we use 

X27 = X? - X7 - Xl4 - 1, 

X64 = X7X14 - X7 + Xl4 + 1, 

X77 = -X? + X?4 + 2X7X14 + 2X7 + Xl4, (82) 

X77' = X? - X? - 2X7X14 - X7 - Xl4 

to derive the following G'2-Ward identities 

= (u{-2dl + 8du + 2d, + U) 

+ v{7dl - 3dud, + 8du - Id, - 7) - 125„) Zg^u, v), (83) 
= {u{ydl-?>dud, + 8du-7d,-7) 

+ viQdl + ml - Qdl - I2dud, - 229„ - 49, + 14) - 249,) (u, t;) (84) 

for the generating function 



Zg2 (u, = j dfircdig) exp (MX7(fi') + ^^Xi4(5')) • (85) 
On the M-axis these equations simplify to 

= {u{-dl + Adu + d, + 7)-6du)ZG,{u,v)U=o, (86) 
= {u{7dl-3d,d, + 8du-7d,-7)-2Ad,)ZG,iu,v)l=o. (87) 

We solve the first equation and its w-derivative for the f-derivatives at f = occurring in 
the second equation and end up with the third order differential equation 

{u'^dl + (14m - 4:u^)dl + (42 - 18m - 19m^)9„ - 14^^ - 56m) Zg^ {u, 0) = 0. (88) 

Again this is probably not related to a generalized hypergeometric series. It may be solved 
in terms of a series expansion 



Zg,{u,0) = 1 + u' + J2^^'' (89) 



n>3 
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provided that the coefficients satisfy the recursion relation (for n > 2) 

9n+i = I , }, , ^^ f 2n(2n + 7)^„ + n(19n + 37)^„_i + 14n(n - 1)^„_2) , (90) 

together with (y^o = fi'i = 0, (72 = 2. These coefficients are related to the triangulations of 
n-gones with inner vertices with valences > 6 [2lj. 

Starting from Zq^{u^ 0) = Zq{u)^ we can now solve the corresponding characteristic problem 
for the first Ward identity (l83ll (analogously to the Ward identity for Spin(5)) by means of 
an expansion 

00 

ZG,{u,v) = Y,~^n{u)V^ (91) 



n=0 



with 



\u) 



uZi{u) = {udl-Audu + Qdu~7u)Zo{ 
2uZM = (l-n)(792 + 89,-7)Z„_2( 

+ {2udl - 8udu + 3{n + 3)du - Uu + 7n - l)Zn-i{u) (92) 

for all n > 2. With the help of (l90ll and ( l92ll . one can reproduce the moments given in the 
appendix (section [B]) recursively. 



5 Results for 5/7(4) 

In this section, we will derive a solution to the Ward identities of the rank 3 group SUi^A) 
in a certain range of the parameter values. For SU{A), the quadratic Casimir of the repre- 
sentation with highest weight = [ni^n2,n^ is given by 

= ^(3?T-i + + ?>n\ + Anin2 + 2nin^ + An2n^) + Srii + 4^2 + 3^3, (93) 

and the fundamental representations 4, 6, 4 with highest weights /xi = [1, 0, 0], /i2 = [0, 1, 0], 
and /is = [0, 0, 1] have Casimirs 

15 , , 

C4 = Q = — and ce = 5. (94) 

The real representation 6 coincides with the vector representation of S'0(6), and 4 is complex 
conjugated to the defining representation 4; the latter two can be identified with the complex 
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fundamental spinor representations of Spin(6). Their tensor products can be decomposed 
according to 

404 = 60 10, 606 = 1© 20 ©15, 4©4 = 6©l0, 

4©4 = 1©15, 4©6 = 4©20', 4©6 = 4©20', (95) 

where again we denoted the representations by their dimensions, 

10 = [2, 0, 0], To = [0, 0, 2], 15 = [1, 0, 1], 20 = [0, 2, 0], 20' = [1, 1, 0], 20' = [0, 1, 1]. 

The representations 15 and 20 are real and 4, 20' are complex conjugate to 4, 20'. From 
((95]), we find 

Xio = xl- X6, Xi5 = X4X4 - 1, X20 = xl- X4X4, X20' = X^Xg - Xi- 
These data enter the Ward identities (l35l l for the generating function 

Zsui^{u, v,w) = j d/ired e">^^+'^^^+"'>^^ (96) 

which is center-symmetric, Zsu{4:) {iu, —iv, —w) = Zsu{4:){u,v,w). They take the form 

= {I5du + ui3dl - 8d^) + 2w{dA - 6d,) + v{dud, - 16)} Zsu{4){u,v,w), (97) 
= {lOd^ + u{dud^ - 6d.,) + 2widl - 2dud, - 4) + v{d,d^ - 69„)} Zsu(4){u, v, w), (98) 
= {159, + v{3dl - 8d^) + 2w{d,d^ - Qdu) + u{dud, - 16)} Zsu{a){u, v, w). (99) 

In order to find an explicit solution for arbitrary products of 4-dimensional representations 
(i.e., on the diagonal with u = v and w = 0), we proceed analogously to section [4Tl We act 
with the operator {udu + vd^ + 9) on the first and last equations and use f l98l) to eliminate 
the w-derivative in the resulting differential equations. At w = 0, we obtain 

{^u'dl + Auvd% + v^d^dl + A^udl + 2bvd^d:) Zsu{A){u, v, 0) 
+ ((135 - 6Auv)du - (IGv^ + A8u^)dy - I60v) Zsu{4)iu,v,0) = 0, 

{3v^dl + ^uvd^dl + n'dld, + ^bvdl + 25u9„9,) Zsc/(4)(m, v, 0) 
+ ((135 - 64mi;)9, - (16^^ + ^8v^)du - 160u) Zsu{^){u, v, 0) = 0. (100) 

Acting with the operator {udu + 2ud^ + 7) on the first equation in fllOOp and evaluating the 
result at M = f = t and w = 0, we obtain the following differential equation for Z(t, t, 0), 

(t^— + 24^2— + (165 - 64t2)t_ + 9(35 - 64^^)- - 960t jZst;(4) = 0, (101) 
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where the derivatives with respect to t are given by expressions analogous to (l55l) . Restrict- 
ing Zsu{4) to the diagonal breaks the Z4 center symmetry down to Z2; we identify (jlOip in 
terms of the Z2-invariant coordinate y = 16t^, 

{y'd'y + 15y'd'y + {QOy' - y')dl + (60y - hy^)dy - ^iy)ZsuiA) = 0, (102) 

as the defining equation f l2Tl ) for the hypergeometric function 

-3/2,5/2 



^5C/(4)(^,^,0) = / rf/ircd(fi')exp (t(tr5( + trsf'f)) = 2i^3 



3,4,5 



(103) 



This result proves the conjecture in [H] (which is based on numerical observations). As a 
by-product, this also leads to the remarkable identity 



^det 



/ In 


4+1 


4+2 


4+3 \ 




In 


4+1 


4+2 




In-l 


4 


4+1 


\ln-3 


4-2 


4-1 


4 / 



(2t) 



2-^3 



3/2,5/2 
3,4,5 



I6r 



(104) 



cf. eq. ( fT6l) for M = t> = t, relating generalized hypergeometric functions and determinants 
of Bessel functions. 



6 On the reduced Haar measure 

In this section, we will describe an alternative approach to Ward identities for the generating 
function Zq based on a factorization of the reduced Haar measure on the maximal Abelian 
torus in G. Tangent vectors to this torus are linear combinations = f pHp of the 
Cartan generators Hp. The reduced Haar measure dfi^ed = PmA dJ'v the maximal abelian 
torus has the product representation [25] 

Pred (e^^'i oc n4sin2 {^-a{HS) = \{ 4sin2 (^-{m,K^)\ (105) 

a>0 ^ ^ mG*+ ^ ^ 

with one factor for every positive root a. The |(dim(G) — rank(G)) positive roots are linear 
combination of the simple roots, 

a = mia(i) + m2a(2) + . . . + "^^^(r), mj G No, (106) 

and the range for m = (mi, . . . , m^)* in (|105l) is chosen in such a way that it parametrizes 
all positive roots. We may take the square root of the density [25] , 

Pred [e'^^) oc |A|2, A = JJ 2isin (-{m,Kip) \ = ^ sign(u;)e™('')(^^), (107) 
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where the sum runs over the Weyl orbit W of the Weyl vector p introduced in (l39l) . Since 
the Wevl orbit of p contains \W\ elements the product representation is preferable for large 
groupsO But it is evident from the sum representation that A changes sign under Weyl 
reflections. 

The density p^-ed oc AA of the reduced Haar measure is a Weyl-invariant function on the 
maximal Abelian torus and hence a function of the fundamental characters. From fllOSp we 
see that it actually is a polynomial of the fundamental characters. In contrast, A is not 
Weyl-invariant and hence cannot be written as function of the characters. 
Less obvious is the observation that A is related to the Jacobian of the transformation 
^ ^ xi'f) from the angular variables to the fundamental characters, 

J{X) = I det ( 1^ ) I oc |A|, such that dp.^a = J^d^'if = J{x)d''x- (108) 



This mapping is one-to-one on the fundamental domain JF of the action of the Weyl group 
on the maximal Abelian torus. This is the closed connected region containing x = in 
which A > 0. We have no proof of the conjecture ( llOSp for all compact simple groups, but 
have checked it for the groups SU{3), Spin(5) and G2 considered in the following sections 
as well as for SU{2) and S't/(4). 

Based on this conjecture we may derive alternative Ward identities from 

/ d^x^{j\x)F{x))=0, (109) 

where use was made of the fact that the Jacobian vanishes on the boundary of the funda- 
mental domain. This leads to the general and simple looking Ward identities 

for any regular function F = F{x) on the fundamental domain. In particular they imply 
the following differential identities for the generating function Zc{u): 

-^(9) + u,J\d)j Zg{u) =0, p = 1, . . . , r. (Ill) 

These should be compared with the geometric Ward identities (l35l ). In reference [17] by 
one of the authors, both types of identities were applied to calculate effective Polyakov loop 



^For example for SU{N) the sum has A^! terms, whereas the product has only iN{N — 1) factors. 
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dynamics of SU{3) Yang-Mills theories on the lattice. The geometric Ward identities are 
usually simpler but not necessarily favored in computer simulations. 

Ultimately, the two systems of linear partial differential equations (I 111 I) and (l35ll must be 
equivalent; but a proof is not straightforward at all. For example, for SU{3) the equations 
(11111) are 4th order differential equation whereas (l35ll are of 2nd order. 
Now we apply the general result (11111) to all simple compact simply-connected groups of 
rank 2. As mentioned above, for these group the conjecture that the density of the reduced 
Haar measure is proportional to the square of the Jacobian J of the transformation ^ X 
can be checked explicitly. The Jacobians for the three groups are computed in the following 
subsections. 



6.1 The group SU{3) 



As in section [4?n [1, 0] denotes the defining representation 3 and [0, 1] its complex conjugate 
3. The reduced Haar measure reads 



red 



1 



6^2" ^ Q^2"v-^-— 3, (112) 

with oc pred ffom (11051) . As a function of the fundamental characters the Weyl-invariant 
and center-symmetric reads 



2r + X3 + Xs-^(9 + X3X3)' 



(113) 



It can be easily verified that its positive square root J indeed coincides with the Jacobian 
of the map (<^i,v92) ^ (X3,X3)- 

The Jacobian vanishes for 



-^X3 




y 



-(9 + 12x + a;^)±2(2x + 3) 



3/2 



= (3{6dl-9d,-d^dl) 



where x = ^{xs) ^ind y = '^{xs) ^ire the real 
and imaginary parts of Xs- The fundamental 
domain inside the triangularly shaped region 
is symmetric under Z3 center transformations 
which rotate X3 by multiples of e^'^*/^. Its 
corners are the values of (9^X3)^X3) at the 
center elements. Here, the identities (11111) 
for the generating function take the form 

u{27 + 4dl + - ISd^d, - dldl))Zsui3), (114) 
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= (3{6d',-9du-d,dl)+v{27 + 4dl + 4dl-18d^d,-dldl)^Zsui3). (115) 
In contrast to the geometric identities (Il8l HOl) these are 4th order differential equations. 



6.2 The group Spin(5) 



With the conventions used in section 14.21 [1,0] = 5 is the defining representation and 
[0, 1] = 4 denotes the spin representation. The reduced Haar measure reads 



with Jacobian J such that 



= (3 + X5 - 2X4)(3 + X5 + 2X4)(4 - 4^5 + xD- 



(116) 



(117) 



The zero locus of the Jacobian is given by 



(-3,0) 








2?/ = ±(a; + 3), y = ±Vx - 1, 

where we abbreviated x = xs and y = Xi- 
The fundamental domain inside the triangu- 
larly shaped region is invariant under the Z2 
center symmetry flipping the sign of Xi- Two 
of the corners, (^5(1), ^4(1)) = (5,4) and 
{x^{z),Xi{z)) = (5,-4), respectively, are lo- 
cated at the values of the characters at the 
center elements. The identities (II 111) for the 
generating function take the form 

u{{3 + - Adl) (4 - 49, + d^) + 3{dud^, - Qdl + lld^, - 209, - 6)) Zspin(5),(118) 

v{{3 + du)^ - Adl) (4 - Adu + d^,) +3{dld, - 8d^, + 229,9, - 79,)) Zspin(5). (119) 




Again, these are 4th order differential equations. 



6.3 The group 



With the conventions of section l473l 7 = [1, 0] denotes the 7-dimensional representation and 
14 = [0, 1] the adjoint representation. The density of the reduced Haar measure 

dfij-ed oc J'^dipid(f2 = Jdx^dxi (120) 
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is a quintic polynomial in Xr ^ cubic polynomial in Xi4, 

= (4x? -X7- 2X7 - 10X7X14 + 7 - 10x14 - Xm) (7 - X? - 2x7 + 4xi4) • (121) 

Since the center of G2 is trivial, this polynomial shows no symmetries at all. Nevertheless, it 
is possible to characterize the fundamental domain for the exceptional group G2 explicitly. 

The zero locus of the Jacobian is given 

by 




y 
y 



(x + l) 



-5(a; + l)±2(a; + 2)3/^ 



where we introduced x = X7 ^ind y = 
XiA- The fundamental domain is the re- 
gion bounded by the three curves defined 
by the above equations. The upper right 
corner is located at the characters of the 
unit element, (x7)Xi4) = (7,14). In this 
case, the identity (11 111) for the generat- 
ing function reads 



(u{Adi - dl - 2du - 109,9, + 7 - 109, - dl) (7 - 92 - 29, + 49,) 
+ 3(3 - 299, + 139;; + 139^ - 389,9, - 21 - 479, + 9^9, - 69,')) Zg^, (122) 
(^;(49^ -dl- 29, - 109,9, + 7 - 109, - dl) (7 - dl - 29, + 49,) 
+ (3 - 299, + 1393 + 1392 - 389,9, - 21 - 479, + 9^9, - 692))^^,. (123) 



Since is quintic we arrive at complicated 5th order linear partial differential equations 
which should be compared with the equivalent geometric Ward identites (l83l IMl) . 

7 Recursion relations for the moments 



The function Zg{u) generates the moments 



"mi ,...,mr 



"/^red Xl ■ ■ ■ Xr 



(124) 



22 



by multiple differentiation at u = 0, see ([9]). Due to center symmetry of the Haar measure 
only center symmetric moments are nonzero, and this selection rule must be respected by 
any recursion relation for the moments. 

One may use the Ward identity for Zg{u) to find such relations. A more direct derivation 
takes advantage of (l33ll with F = xT^ ■ ■ ■ xT'' ■ finds 

q 

+ J2 Cp,c,m, I rf/i,ed XA(x)xr ■ ■ ■ XT~' ■■■xT- (125) 

Clearly, the complexity of these relations increases with the degree of the polynomials xx 
in the last sum. Alternatively, we could apply eq. (11 101) to the same function F, with the 
result 

= j (i/i.ed (3 ^^^XT ■ ■ ■ Xr + 2m,j\x)xT ■ ■ ■ X7-' ■ ■ ■ ) • (126) 

These relations are based on the conjecture (jlOSp . in contrast to the 'geometric recursion 
relations' in (I125p . 

For the group SU(3): For this group the recursion relations p25l) take the form 

(8 + 2m + n)t„,^x,n - 6mtm-i,n+i - 9ntm,n-i = 0, (127) 
(8 + 2ra + m)tm,n+i - Qntm+i,n-i - 9mtm-l,n = 0. (128) 

Since the moments are symmetric the two identities are equivalent. These 'geometric iden- 
tities' are much simpler than the 'non-geometric' relations (I126p . which for p = 1 lead 
to 

(18 + 4:m)tn,m+2 - (3 + m)tn+2,m+l - (27 + 18m)t„+i,„ + 27mtn,m-l + 4:mtn+3,m-l = 0. 

By symmetry of the coefficients tmn, the relation for p = 2 is again equivalent to this 
recursion formula. The difference of both leads to the simpler formula 

+m+3,m 

+ {Ak - 6)t3k+m,m+3 
+ k (27t3fc+m,m ~ 18t3fc+m+l,m+l ~ ^3A;+m+2,m+2) • (129) 

All recursion relations are compatible with center symmetry which implies = unless 
m = n mod 3. With the moments 

2(3m)! 6(3m+l)! 
" m!(m + l)!(m + 2)! ' " m!(m + l)!(m + 3)! ^ ' 
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one can compute all tmn with the recursion relation (I127p . For example, for m = n, one 
finds 

^ f2k\ fm\ ^ 3k^ + 2k + l- 2km -m 
t^m - ( ^ J (^^ J (A; + l)2(A; + 2)(m-fc + l)- ^^^^^ 



k=0 



The moments tmn for small m and n are given in the appendix. 

For the group Spin(5): For this group the geometric recursion relations (I125p read 

(8 + 2m + n)tm+i,n + (4m - 5n)tm,n - 6mt„„i,„ - 4:mtm-i,n+2 = 0, (132) 
(5 + m + n)tm,n+i - 5mtm-l,n+l - 2ntm+i,n~i - Qntm,n-i = 0. (133) 

These relations are compatible with center symmetry which implies that tmn = for odd 
n. With the help of the first or second recursion relations, one can determine any tmn given 
tmfl- For the first relation we also need the to,n which are just the coefficients in the series 
expansion of Z{0,v) in (iTOl l. The moments tmn for small m and n can be found in the 
appendix. 

For the group G2: For this exceptional group the recursion relation p25p are more 
involved, 

= 2mt„_i_„+i - (12 + 2m + 3n)tm+i,n + (8m - 7n)tm,n 

+ lAmtm-l,n + 7ntm+2,n-l + 8ntm+l,n~l 7Titm,n—li 

(134) 

= (24 + 3m + 6n)tm,n+i + 7mtm-i,n+i - (7m - 12n)t„+i^„ - (8m - An)tm,n 

+ 7mtm-l,n — Qntm+3,n-l — '^0ntm+2,n-l + 22'ntm-\-l,n—l 14?^tm,n— !• 

(135) 

For example, one can calculate all tmn from the tm,o and tm,i- The former ones are just the 
coefficients gm in (l90ll . The moments tmn for small m and n are given in the appendix. 



8 Conclusions 

In this paper, we have derived two kinds of Ward identities for the generating functions 
for integrals over arbitrary polynomials of fundamental characters. One is a consequence 
of the fact that a left-derivative of any function on the Lie group integrates to zero with 
the full Haar measure. For a convenient choice of this function, this left-derivative is a 
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class function so that the vanishing of the integral reduces to an identity on the maximal 
Abelian torus. The other Ward identity derives from an integral of a total derivative of an 
arbitrary class function over a certain domain. If one chooses this domain as the region 
where the Jacobian of the change of variables from the angles of the maximal Abelian torus 
to the fundamental characters is non-negative, one can split powers of the Jacobian from 
the arbitrary class function so that the result vanishes on the boundary of this fundamental 
domain. This, however, leads to ultimately more complicated differential equations for the 
generating functions than the first, geometric, approach. Both furnish generalizations and 
structural clarifications of identities used in the case of SU (3) in an earlier publication [T7] 
by one of the authors. In this paper, they have been applied to all simple compact simply 
connected Lie groups of rank two. Furthermore, they have been used to prove several 
conjectures in the literature concerning explicit solutions for SU (3) and SU (4); beyond that, 
we derived recursion relations determining all integrals over polynomials of fundamental 
characters for the above groups. 

The derivation of the second kind of Ward identities is based on a conjecture concerning 
the factorization of the reduced Haar measure density pred into the square of this Jacobian. 
This conjecture has been checked explicitly for several cases under consideration, but so far 
it lacks a general proof. It might be interesting to clarify this issue from a group theoretical 
point of view. 

Another obvious open question is the equivalence of our two approaches. Since both encode 
information determining the same generating functions it ultimately should be possible 
to derive one from the other. Perhaps with a deeper insight into the group theoretical 
connection between the approaches it might be possible to obtain even simpler identities 
which might be of use, e.g., in lattice gauge theories or random matrix models. From a 
mathematical point of view they answer the question how many invariants there are in a 
given tensor product of fundamental representations. 

We found it surprising to note that the generating function for powers of the sum of char- 
acters of only the defining and its complex conjugate representation for SU{2) and S't/(4) 
can be expressed in terms of appropriate generalized hypergeometric functions. One might 
speculate that this generalizes to higher rank SU{2n) as well. Apart from that, this fact 
for SU{4) leads to a hitherto unknown relation (11041) between generalized hypergeometric 
functions and Bessel functions. Finally, one might note the remarkable fact that all Ward 
identities (of the first kind), at least for the Lie groups of rank 2, reduce to parabolic dif- 
ferential equations on the locus where their families of characteristics coincide. This allows 
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for a recursive integration starting from a solution on this locus. 
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A Generating function for SU{3) 

In this section, we derive an analytic solution to eq. (I48ll . By invariance of the Haar measure, 
dfi{g) = dfi{g~^), we expect Zsu(3) to be a symmetric function of u and v. Thus, our 
objective is to solve (l48l ) and (l49l l with the symmetry Zsu{3){u,v) = Zsu{:i)iv,u) and initial 
condition Zsu{S){.^i^) = 1- From fl23l) . we know that 



Multiplying (1481 ) by u and ( l49l) by f , the difference of the resulting equations reads 



In principle, it is possible to find a recursive solution analogously to (TfTl) and (lOTl) starting 
from a solution on the locus where the families of characteristics of the differential equation 
coincide (here on the u- and the v-axis), where the differential system turns out to be 
parabolic. However, in the case of 5'f/(3) we are able to give a closed expression of the 
solution which allows for an easier computation of the moments ([9l): 

Center symmetry (ISOll as well as the symmetry of Zsu{2,) in its two arguments suggests to 
introduce new coordinates x = uv and y = + v^. As a function of these variables, Zsu{3) 
has to satisfy 




{u^dl - v'^dl + 4 {udu - vd,) - ^{u^d, - v^du))Zsu{3) = 0. 



(136) 



(3?/9j + Qdy + 2xd^dy - d:,)Zsu{3) = 0, 



(137) 



which can easily be solved by a power series expansion, 



oo 




(138) 



m,n=0 



with recursion relation 



{3n + 3 + 2m)na 



mn 



(m + l)am+l,n-l- 



(139) 
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The solution to (I138P and (I139P with the condition 



oo 2 °° 

n=0 ^ ^ ^ ' ^ n=0 



is given by 



Zsui3){u,v)= M H M^ + t;^ . 141 

^-^ n! m + n + 1 ! m + n + 2 ! V ml 

m,n=0 ' 

This solution also satisfies the original diff'erential equations (Il8l) and (ll9l) . It encodes 
information about how many invariants there are in the tensor product 3®^ ® of 5'[/(3) 
representations. General results for integrals over S'f/(3) matrix elements can be found 
in [121. 



B Tables of moments for the rank 2 groups 

With the recursion relations ffT27l) . ffT28D . ffT32D . (fT33D . (fT3iD . and (fT35ll . we determined the 
moments 

tm„= / C?yUrcd(^)X[5',o](^)xfo,l](^), (142) 

Jg 

for the groups with rank 2 for small m and n. Here X[i,o] and X[o,i] are the characters of the 
fundamental representations with highest weights \i\ = [1,0] and fj,2 = [0, 1]. 
For SU(3) the lowest moments of xTx^ are 



m\n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





1 








1 








5 








42 





1 





1 








3 








21 








210 


2 








2 








11 








98 








3 


1 








6 








47 








498 





4 





3 








23 








225 








2709 


5 








11 








103 








1173 








6 


5 








47 








513 








6529 





7 





21 








225 








2761 








38265 


8 








98 








1173 








15767 








9 


42 








498 








6529 








94359 





10 





210 








2709 








38265 








586590 
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For Spin(5) the lowest moments of x^x" are 



m\n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 





1 





1 





3 





14 





84 





594 


1 








1 





5 





30 





210 





1650 


2 


1 





2 





11 





75 





580 





4917 


3 








4 





27 





205 





1714 





15435 


4 


3 





10 





73 





600 





5338 





50506 


5 


1 





26 





211 





1852 





17342 





171022 


6 


15 





75 





645 





5970 





58350 





596085 


7 


15 





225 





2061 





19950 





202230 





2129719 


8 


105 





715 





6837 





68730 





718928 





7774600 


9 


190 





2347 





23403 





243050 





2612796 





28922112 


10 


945 





7990 





82301 





879204 





9681144 





109404729 



Both tables nicely display the constraints from center symmetry. 
For G2 the lowest moments of xTXu 



m\n 





1 


2 


3 


4 


5 


6 


7 


8 





1 





1 


1 


5 


16 


80 


436 


2786 


1 











1 


6 


40 


260 


1785 


12852 


2 


1 


1 


3 


10 


45 


236 


1421 


9444 


67852 


3 


1 


2 


7 


32 


170 


1016 


6637 


46656 


348553 


4 


4 


9 


33 


151 


817 


4984 


33357 


240181 


1835171 


5 


10 


30 


126 


641 


3728 


23986 


167080 


1241285 


9727650 


6 


35 


120 


545 


2932 


17827 


118945 


854135 


6511050 


52159514 


7 


120 


476 


2359 


13517 


86171 


596686 


4415055 


34500369 


282217558 


8 


455 


2002 


10626 


64078 


425194 


3041241 


23115050 


184754906 


1540766892 
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